A 233 km Circumference Tunnel for e+e , pp, and p^p Colliders 
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ABSTRACT 



In 2001 a cost analysis survey was conducted to build a 233km circumference tunnel in 
northern Illinois in which to build a Very Large Hadron Collider. Ten years later I have 
reexamined the proposal, taking into consideration the technological advancements in all 
the aspects of construction cost analysis. I outline the implementations of e+e", pp, and 
//"^/i" collider rings in the tunnel using 21st century technology. The e'^e" collider employs a 
Crab Waist Crossing, ultra low emittance damped bunches, 12 GV of superconducting RF, 
and 0.026 Tesla low coercivity grain oriented silicon steel/concrete dipoles. The pp collider 
uses the high intensity Fermilab p source, exploits high cross sections for pp production of 
high mass states, and uses 2 Tesla ultra low carbon steel/YBCO superconductor magnets 
run with liquid neon. The /x^/x" ring ramps the pp magnets at 8 Hz every 0.4 seconds, 
uses 250 GV of SRF, and mitigates neutrino radiation with a phase shifting roller coaster 
FODO lattice. Such coUiders could provide the means necessary for Chicago, and more 
importantly America, to stay relevant and competitive in the international marketplace of 
particle physics. 
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Introduction 



For many years the University of Mississippi has collaborated on experiments in the 
Fermilab fixed target [H [2] and collider [21 S] programs. In this thesis, we will explore the 
design of a large, 233 km circumference multipurpose ring for energy frontier colliders, uti- 
lizing more cost effective machines and methods. The colliders use several new technologies 
developed and refined in the last decade. 

Advances in the automation of some, but not all, aspects of tunnel boring machines and 
wall reinforcement are explored. Improvements to the tunnel's structural integrity will also 
be discussed. 

For the e~^e~ collider very low emittance bunches are combined with a relatively large 
crossing angle and a short vertical focal length at the interaction point for high luminosity 
by utilizing new methods of focusing collision regions implementing the Crab Waist crossing 
design. Also, we can integrate the use of low coercivity, grain oriented silicon steel for dipoles. 

A pp hadron collider is chosen because the cross section for many high mass objects 
is much larger than for pp collisions. Single bore iron dipoles permit high temperature 
superconductor ribbons running at liquid neon temperatures. Advances have been made in 
p production at Fermilab. A detector with more physics and fewer background interactions 
per crossing requires fewer channels and will consequently cost less. 

The same dipoles used for pp can also be ramped to accelerate muons to almost 20 TeV 
with less RF than the proposed International Linear Collider. The ramping rate is slow 
enough to permit superconductors. A simple method which solves the neutrino radiation 
problem will be worked out. A 40 TeV center of mass muon collider with an effective energy 
reach 30 times higher than the CERN LHC is envisioned to do physics [S] |6] at Fermilab. 
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Chicago is one of the great cities in the world and it needs to be home to one of the 
greatest scientific instruments in the world to maintain its position on a competitive planet. 
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1 Boring a 233 km Long Tunnel: 
A Construction Cost Analysis 



1.1 Drilling 

Essentially in the world of major underground drilling there are three major techniques 
commonly employed. These methods include tunneling machinery that implements direc- 
tional drilling bits, microtunnel boring machines (MTBM), and tunnel boring machines 
(TBM). The most basic and commonly used of the three is directional drilling. Applica- 
tions of directional drilling are vast, but seen mostly in oil well preparation. The large 
advantages of directional drilling techniques include a relatively low cost of operation (as it 
is unmanned), a multi directional drilling path, and the most easily excavated milling bi- 
product (a fine silt/water mixture referred to as "muck"). Also the need to change out the 
tunneling bit head has greatly reduced as synthetic diamond tipped and diamond ceramic 
bit heads have become increasingly available. The overwhelming disadvantage of directional 
drilling however, comes in the form of a very limited drill head diameter, topping out at 
approximately 900mm; however, this might still be large enough for access shafts. For the 
purposes of this paper, we will focus most directly on tunnel boring machines (TBM) and 
microtunnel boring machines (MTBM). 

Like directional drilling, microtunnel boring machines are most commonly operated 
through remote control (thus largely unmanned). While this is also a great advantage to the 
cost of operation and efficiency of a tunneling operation, MTBMs also share the limitations 
of directional drilling in the fact that they too have dimensional boring size restraints. While 
MTBMs are definitely larger than directional drilling methods, topping out at approximately 
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1.5 meters ( 5 ft), there would still not be enough space for manned transportation within 
the underground collider ring. Also, most conventionally MTBMs incorporate a pipe jacking 
mechanism that essentially drives a preassembled pipe into the tunnel to reinforce the walls 
of the tunnel. Problems occur in this model as after a few kilometers of pipe jacking, fric- 
tion does not allow for further pipe lengths to be installed. However, this tunnehng method 
should not be completely discounted, as it might be useful for installing larger elevator shafts. 

Tunnel boring machine methods are by far the most useful methods introduced through- 
out this section. TBMs are the most costly method of underground tunnehng excavation, 
but are also the most sizable. Often referred to as a "mole" , TBMs have large circular faces 
that rotate and dismantle any solid terra in its path. By utilizing a pressurized mechanism 
that breaks away the rock face while allowing the large chunks of rock to fall through open- 
ings in its rotating face, TBMs are extremely efficient in burrowing through and disposing 
of solid surfaces. One obvious advantage of breaking apart the tunneling face into larger 
chunks of particulate is the vast amount of energy saved in not grinding and refining the 
particulate down into a finer milled product. In acquiring larger chunks however comes the 
problem of disposing of a larger milling byproduct. These larger pieces of rock typically have 
to be disposed of by an integrated conveyor belt system that moves the rock away as it falls 
through the grinding mechanism face. Although there are many methods of reinforcing the 
tunnel while the drilling continues throughout its directed path, most methods incorporate 
some combination of drilling pilot holes in the walls of the tunnel and placing reinforcement 
bolts with wire mesh, solid walls, concrete matrices, etc. to eliminate the possibility of wall 
collapse. The biggest advantage of utihzing TBM methods comes in the sheer size of the 
tunnel TBMs can accomplish. Though most TBM faces are 3-5 meters (~10-16 feet) in 
diameter, current TBM manufacturing companies are pushing nearly 16-meter (~52 feet) 
production models. 
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1.2 Past Cost Estimation 

In 2001 an initial cost estimation was projected for the construction of a very large hadron 
collider tunnel (VLHC) to be constructed in Northern Illinois for the Fermi National Accel- 
erator Laboratory. This report was prepared by CNA Consulting Engineers and prepared an 
array of building initiatives and a very comprehensive cost analysis for several scenarios in 
which the collider tract could be constructed [3 [8]. The final cost projection put the entire 
project in the $2.5-2.8 billion dollar range, depending on tract orientation and location. All 
of the intrinsic costs will be further explored in the following paragraphs. 

One of the most fundamental parts of estimating the cost of tunneling underground comes 
in a simple equation that determines the quality (Q) of the material being drilled. Quality 
is most directly correlated to the hardness and softness of the rock being drilled. The harder 
the rock, the higher the Q factor, thus the better the quality. Conversely, as rock is softer, 
its Q goes down, and its relative quality is considered bad. Typically the range of Q is set 
somewhere between .2 and 50 and is determined by the following equation: 

Q = ^ — ^ -F, — (1) 

Jfi J a SRF 

Where from Eq{T] , RQD = Rock Quality Designation , J„ = Joint set number, = Joint 
roughness number, Ja = Joint alteration number, Jy^ = Joint water reduction number, and 
SRF = Stress Reduction Factor, respectively. Although all these variables are representative 
of physical characteristics of the rock's internal physical composition, the overall Q factor 
really reflects the structural rigidity and integrity of the rock to resist collapsing on itself. The 
most cost efficient option for building this collider tunnel came in the form of an inclined track 
built in the northern portion of the state. The cost effectiveness of this specific projection 
came in the limited range of Q quality factor (< 33) and the effects of inclining the track itself, 
so that major bedrock drilling could be put to a minimum. The projected range of Q quality 
for this tunnel to be pushed through the northern incline ranged from ~. 33-33. Overall, this 
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ranging factor essentially represents the number of reinforcement measures needed per unit 
of tunnel. For example, in the low end of the spectrum (0.33), reinforcement rock bolts would 
need to be drilled every 1.5 m and the thickness of the shotcrete wall covering the tunnel 
surface would need to be around 0.22 m thick. Conversely, on the other end of the Q factor 
spectrum, the most rigid well (~33) could be composed of rockbolt spacings approaching 
2.6 m and shotcrete thicknesses of around 0.07 m. So it is easy to see that the integrity of 
the wall itself plays a large role the construction cost of the entire project. 

It is also worth noting the cost per square foot of cross section of each projected tunnel 
construction in past projections. When considering the 12 ft diameter tunnel at a circumfer- 
ence of 233 km, roughly 2,446,500 m^ of tunneling space would be created. If we extended 
the interior tunnel diameter to 16 ft, at the same circumference, roughly 4,350,000 m^ of tun- 
neling space would be created. So by extending the interior diameter of the tunnel from 12 
ft to 16 ft you would essentially double the volume of the tunnel; however, the cost of the 16 
ft tunnel is not double the 12 ft. In fact the 12 ft construction cost came in at approximately 
$2.42 billion, while the 16 ft came in at $2.71 billion. A difference of approximately 10% 
cost difference compared to a 100% bigger tunnel. Roughly 30% of the total construction 
costs come from contingency costs. A large part of this is the uncertainty of the actual costs 
of manning the mole operation. 
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Figure 1: Cost Analysis of Various Tunnel Diameters and Orientations [7j 
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1.3 Future Projections 

Throughout examining the past projections of construction costs, many improvements in 
drilhng technologies and cost effectives have arisen. In the past decade great advancements in 
unmanned drilhng possibilities have become evident, as well as the overall increases in diam- 
eters and efficiencies of directional drilling techniques, TBMs, and MTBMs. Some of these 
advancements will be addressed in the following section as well as some new construction 
techniques that could serve for another building model altogether. 

One such construction tactic could be the implementation of periodic reinforcement bulk- 
heads throughout the shotcrete lined tunnel. As it stands in past projections, shotcrete 
thicknesses ranged from 0.07— 7-0.22 m in the lining of the completed tunnel. This thickness 
differential is essentially varying by a factor of 3 and if it could be reduced to a thickness 
of 0.07 throughout the entire tunnel, roughly 200,000 m^ of shotcrete could be saved. This 
might be accomplished by installing period steel bulkheads to reinforce the structural in- 
tegrity of the tunnel, much like the construction of modern day submarines. The structural 
rigidity of a thin walled cylinder to resist external pressure increases by enclosing the ends. 
The resistance to collapse is derived from the equation for collapsing pressure, [9], 

= X E X psi (2) 

Where from Eq|2], K is a numerical coefficient that contains the stress variable, E is the 
elasticity modulus of the material, D = outside-shell diameter, and t = the thickness of the 
shell. The modulus of elasticity of steel is ~200 GPa compared to ~20 GPa for concrete. 
Adding steel bulkheads periodically, as in a submarine, to the thin shotcrete lined cylindrical 
tunnel would save a lot of unnecessary shotcrete and add a lot more structural integrity to 
tunnel. Also, the idea of creating these bulkheads using concrete might seem more feasible, 
but if you consider a 0.5 m thick steel bulkhead versus a 0.25 m thick concrete bulkhead, the 
steel bulkhead would be 80 times more resistant to buckling. Obviously this would mean far 
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fewer bulkheads, thus construction costs would be further reduced. 

The advancements in TBM automation are also very interesting when considering the 
construction cost of the VLHC in Northern Illinois. Of the $2,418 billion dollar construction 
budget for the 12' diameter northern inclined track, roughly 30% came from a contingency 
budget. A large percentage of this portion of the budget arises from the sheer cost of manned 
operation and liability of the TBM. With great advancements in TBM automation, many 
TBMs are now operated by relatively few operators. When considering the action of a 
TBM, four distinct operating actions must be considered: 1) The guidance of the TBM, 
2) The setting of stabilizer supports on the tunnel walls for which the TBM establishes 
forward pressure, 3) The setting of wall reinforcement measures, and 4) The removal of rock 
and particulate and their disposal. Through the great advancements in robotic controlling 
technology over the past decade, guidance mechanisms can be mostly operated from the 
surface, with the occasional subsurface manned intervention. The setting of stabilizer legs 
for the forward pressure of the TBM are now completely automated by the TBM itself. 
Essentially, they're just huge retractable legs at the rear of the TBM that continuously retract 
and reset as the TBM moves forward. The setting of tunnel wall reinforcement has become 
almost completely mechanized and automated as well. TBMs can now drill and place rock 
bolts periodically as needed and secure reinforcement mesh. The only reinforcement tactic 
that still requires manned operation is the shotcrete application, and it is being rapidly 
addressed with new automated machinery that might function as a shotcrete robot, that 
would also be unmanned. While rock and particulate disposal still need some manned 
assistance, long conveyor tracts may be attached to the rear of the TBM to pull rock back 
to a more secure section within the tunnel. 

Further investigation could be made into the construction cost analysis of the VLHC 
proposed almost a decade ago. With the continued advancement in TBM automation, maybe 
soon we could see a machine that can bore 8 ft diameter tunnels (a very usable size) without 
much human interference or manned guidance at all. Other methods might be found to 
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minimize costs on a 8 ft diameter tunnel's construction, as an 8 ft tunnel has roughly one- 
fourth the volume of a 16 ft tunnel. 
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2 Electron-Positron 



2.1 Theory 

Electron-positron collision can result in the annihilation of each particle and the release 
of a pair of gamma rays, Eq|3l 

e+ + ^ 7 + 7 (3) 

However if both particles have appreciable energies that are at or beyond the mass of the 
carriers of the weak force, W and Z bosons, then the weak force becomes comparable with 
that of electromagnetism and boson production is possible. Although currently the heaviest 
particles realized in this annihilation are — W~ pairs, the heaviest single particle is the 
Z boson. This highly energetic annihilation is the driving force behind the proposal for the 
International Linear Collider, ILC, as it is likely the source of the much sought after Higgs 
boson, {H), FigJl 




Figure 2: Feynman Diagram of Higgs Boson Production. 

2.2 Past Research 

One of the most relevant research endeavors to date was the Large Electron-Positron Col- 
lider, LEP, conducted in 4 experiments at CERN from the early 90's to the early 2000's [lOj . 
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LEP allowed for the production of the Z boson, as well as the pair production of — W~ 
bosons. In 2002, Tanaji Sen and Jim Norem introduced a proposal that would construct a 
Very Large Lepton Collider, VLLC, inside the proposed 233 km collider ring space built in 
conjunction with Fermilab [H]. The parameters surrounding the proposed collider are very 
interesting and the theoretical explanations for this construction will serve as the basis for 
this section. 

Examining the luminosity, L, of a specific Electron-Positron bunch collision, EqJH 



1 



7/ (4) 



2ere f3; 

Where C,y in the vertical beam-beam tune shift, /3* is the focal length in the y plane, 7 
is the relativistic scaling factor (that essentially scales as the Collision Energy), and / is the 
total beam current in both beams. Now if you consider the parameters examined at LEP 
in 1999 and the parameters proposed for the VLLC in Table [H the increase in luminosity 
is easily seen. Essentially you get an increase in the vertical beam-beam tune shift, ^y, of a 
factor of 2.28, a decrease in the focal length /3* of 5, an increase in the relativistic 7 of 2.05, 
and an increase of the total beam current in a both beams, /, of about 3.95. The adjusted 
parameters give you an overall increase in luminosity of 90 times the original luminosity of 
LEP, EqE] 

LvLLC = ^ " i92.31)LLEP = 8.8 x 10=^^ cm-^s-^ (5) 

Although this is a great increase in the overall luminosity of the Electron-Positron colli- 
sions compared at LEP, we can do better and these improvements will be the focus of our 
next section. 

2.3 Future Projections 

One of the most prominent and cutting edge techniques to improve electron-positron 
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collisions comes in the form of creating a Crabbed Waist (CW) crossing at the intended 
coUision site, FigJSl [121 [131 HI]- Instead of having the two beams collide head on, they are 
angled so that the collision only occurs in the section where they overlap. Now the angling 
of beam collision alone would produce a great decrease in luminosity, as it would decrease 
the overlapping collision area. Actually, it would go down exactly by a factor correlating to 
the angle at which the beam is tilted, or more commonly known as the Piwinkski angle, 0, 
where is 

= — (6) 




Figure 3: Crabbed Waist Crossing Scheme [12] 



However, a few more parameters could be tweaked to compensate for the luminosity loss 
and even improve it overall. One of the most important consequences would be the ability 
to lower the vertical focal length at the collision point, /?*, to fit the decreased overlapping 
area. Before we explore any changing of parameters in the luminosity equation though, let's 
explore a more direct approach approximating the luminosity parameters put forth by Sen 
and Norem [TT] . 

Well taking the equation for luminosity to be EqJTl 

12 



The equation for tune shift, ^y, goes as EqJHl 
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y 



(2vr7cry(a^. + ay)) 
Then we can combine EqJ7]and EqJHlto get Eq|9l 

, 7S.iVi/o(l + ?J ,„> 
^ Pi^f)— 

The next step after resolving these equations is to simply put forth EqlTU] [15], 



^ (2.167 X 1034) xE(GeV)x 1(A) x^y 

But where did this constant come from? For the sake of sanity we'll calculate this constant 
directly. For a 1 GeV electron, 7 = ^-^^ = 1957. Taking the classical radius of an electron 
to be Te = 2.818 x 10^^^ cm and the electron charge to be 1.602 x 10^^^ C, then we can 
calculate the constant to be EqJTTl 



^ ^^^^ 2.167x10^4 (11) 



(2 X re X Qe) (2 X 2.818 x lO-i^ x 1.602 x 10-i9) 
plugging this into the parameters outlined in Sen and Norem, Table [H we get Eq 



^ (2.167 X 1034) X (200 GeV) x (0.0114A) x 0.18 ^ ^ ^ .3 _^ _^ 

L = — ^ ^ = 8.8 X 10^3 cm ''s (12) 

(1 cm) 

This result matches the result in Sen and Norem, Table [1], almost exactly [H] . So now 
we can be confident to plug a few more recent parameter values established by more recent 
research endeavors to further our luminosity approximations. One such parameter comes 
from International Linear Collider Reference Design Report, in which a focal length of /3* 
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=0.06 cm is achieved [16]. Plugging this into our current parameters for EqJTOl we get EqJT3] 



^ ^ (2^167 X 10^^) X (200 GeV) X (0^0114A) x 0^18 ^ ^ ^^3. (13) 

(0.06 cm) 

This is astonishing, as luminosity has gone up by a factor of 17; however, to stay in hne 
with current ILC specs, we should really explore the implications of going up to a 250 GeV 
collision energy. One consequence of going up to a higher energy is an increased allowable 
tune shift {E : 200 — t- 250 GeV, corresponds to : 0.18 — )■ 0.23). However, we must consider 
the 100 Megawatt synchrotron radiation limit, so then the current, I {A), scales accordingly, 
(||2)4 = 2.44. The current then goes down by a factor of 2.44, {1(A) : 0.0114 0.00467). 
Plugging in all these new parameters in EqJTO| we get EqJH] 



(2.167 X 10^4) X (250 GeV) x (0.00467A) x 0.23 ^ ^ o i 

^ ^ ^ ^ ^ = 9.7 X 10^^ cm-^s"^ (14) 



(0.06 cm) 

Which is still about 10 times higher than the luminosity from Eqj29]and 5 times higher 
than that of the ILC. Another parameter that needs to be considered are the emittances 
that correspond to the crab waist scheme and their respective collision energies. If we apply 
the crab waist scheme with half crossing angle = 17 mrad, assuming a bunch length = 
6.67 mm, and an overlap length I = .5 mm = then = 4.25 fim. If we know the focal 
length in the horizontal plane to be /3* = 20 000 fim, then the emittance in the horizontal 
plane goes as 



-f = 0.9 nm. (15) 

(3* 



Then calculating the 6N2 term. 
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2N2(T. 

6N2 = 



6N0 



2 X (4.85 X 10^^) X 4.25 urn 



(.0017 X 6.67 mm) 
SN2 = 3.63 X 10^° (16) 



Now we can calculate the collision frequency, /o, for 250 GeV, 



{Number of Bunches / Beam) x c 

Jo = 



(Circumference of Collider) 
(46) X (3 X 10^ ^) 



(233 km) 

/o = 5.92 X 10^ Hz (17) 
Now we can calculate the beam size in the vertical plane, Cy, for 250 GeV, 

msN^fo 

(4.85 X 10i^)(3.63 x 10i°)(5.92 x 10^ Hz) 
~ 47r(4.25 X 10"^ cm)(1.5 x lO^^cm-^s^i) 
ay = .0201 fim (18) 

Just as in EqUS], the emittances in the vertical place, ey, goes as EqJTOl 

crl (.0201 um)2 ^ ^ 

e --^ = i- t^-J- = 0.00067 nm. 19 

^ P; (600 /xm) ^ ' 

We calculated the a^^, ay, e^, and ey for the 200 GeV energy collision as well, utilizing 
the exact same procedure as above. Those values can also be found in Table [H 

One factor that must be taken into consideration is whether or not this emittance is 
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even possible. This is easily done by looking at our emittance projections compared to the 
current standards at the ILC. The current values achievable at ILC for horizontal and vertical 
emittances are = 12 mm-mrad and = .08 mm-mrad ^XQ\. Converting these values using 
our proposed gamma factor for 250 GeV, 7, 



^ ^ EcoiHsion ^ 250 GeV ^ ^^^^^^ ^^q) 
rrip- 511 keV 



Then converting mm-mrad — )■ nm, 



12 X 10-6 m ^^^^^ 

= — = = 0.0245 nm 

7 489000 

e'i 0.08 X 10-6 m , , 

= ^ = = 0.00016 nm 21 

^ 7 489000 ^ ' 

It is easy then to see by Eq|2T]that the values we proposed are both larger than the values 
put forth by the current ILC parameter proposal [e^ ^ 0245° nm ~ ^0, °-°°°67 nm _ 

This is very inspiring because these projections are actually realistic when compared to the 
values from ILC that have already been accomplished with ultra low emittance damping ring 

HZ]. 

When considering that the dipoles for this 233 km circumference ring have a magnetic 
field 4 times lower than used at the CERN LEP machine, to maintain good field quality, a 
soft magnetic material is needed [iHl [iHl |20l [21] . Grain oriented silicon steel is chosen for 
the dipoles because its coercivity is an order of magnitude lower than low carbon steel [22] . 
Horizontal bands are added to the top and bottom of C shaped laminations to permit high 
permeability path in the entire flux return circuit. Putting concrete in between laminations 
would provide space for the two bands. An example of this design is currently being tested 
and calibrated at the University of Mississippi, FigJH |23] . 

The beam-beam bremsstrahlung cross section becomes a fair fraction of a barn, leading 
to 10^^ ■ 10-24 = 10" interactions per second. With this high of an interaction frequency, the 
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beam must be replaced every few minutes. The horizontal tune shifts, K2 strength of Crab 
Waist sextapoles, and maximum beam size in interaction point focusing optics are calculated 
as follows [T2|[2l]: 



29h (3yl3*y 




^^1 1 / 0.02 m 

K2 -- 



2(0.017) (250000 m)(0.0006 m) V 40000 m 
K2 = 0.0014 (22) 

Which appears to be achievable, for a sextapole. Now on to calculating the maximum 
beam sizes in the IP quadrupoles to see if they fit , Eqj23| 



(Tx = \/e^ = VO.9 X 10-9 ■ 40000 = 6 mm 



ay = v^e^ = VO.00067 x lO'^ ■ 250000 = 0.41 mm (23) 
Finally we can calculate the tune shifts in each plane [15], Eqj2i 



7r7 (7^6'^ 



= ^-5^ = 0.027 



^ !£^^ = 0.23 (24) 

Note that luminosity, RF power, and bunch length can be traded off. A longer bunch 
length allows a lower injection energy by minimizing the transverse mode coupling instabil- 
ity [H] . The power for the superconducting RF cavities might come from klystrons or from 
phase and frequency locked magnetrons [25], [26] . 
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Figure 4: 1.8 Tesla, 400Hz Grain Oriented Silicon Steel Dipole Magnet at the University of 
Mississippi [23] 



Table 1: e^e Collider Parameters 



Parameter 


LEP 


VLLC 


CrabWaist2QQ 


CrabWaist25o 


Circumference (m) 


26 658.9 


233 000.0 


233 000.0 


233 000.0 


/S:, [i; (cm) 


150, 5 


100, 1 


2, .06 


2, .06 


Luminosity (cm ^ sec 


9.73 X 10^1 


8.8 X 10^3 


1.5x 10=^^ 


9.7x 10^^ 


Energy (GeV) 


97.8 


200.0 


200.0 


250.0 


7 


191 000 


391 000 


391 000 


489 000 


Emittances e^, (nm) 


21.1, 0.220 


3.09, 0.031 


.9, .0017 


.9, .00067 


rms beam size IP a*, a* (/um) 
Bunch intensity/I (/mA) 


178.0, 3.30 


55.63, 0.56 


4.25, 0.0321 


4.25, 0.0201 


4.01x10" 
0.720 


4.85x10" 
0.1 


4.85x10" 
0.1 


4.85x10" 
0.04 


Number of bunches per beam 


4 


114 


114 


46 


Total beam current (mA) 


5.76 


22.8 


22.8 


9.34 


Beam-beam tune shift 


0.043, 0.079 


0.18, 0.18 


0.034, 0.18 


0.027, 0.23 


Dipole field (T) 


0.110 


0.0208 


0.0208 


0.0260 


E loss / / turn (GeV) 


2.67 


4.42 


4.42 


10.8 


Bunch length (mm) 


11.0 


6.67 


6.67 


6.67 


Revolution frequency (kHz) 


11.245 


1.287 


1.287 


1.287 


Synch rad pwr (b.b.) (MW) 


14.5 


100.7 


100.7 


100.7 
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3 Proton - Antiproton 



3.1 Theory 

One of the most abundant sources of naturally occurring antiproton production occurs 
when cosmic ray protons interact with the interstellar medium. This example will serve as 
a good model for the creation of antiprotons. Essentially you have a very energetic proton 
traveling through space that collides with any arbitrary nucleus in the interstellar medium, 
and the following reaction occurs 

p + A-)- P + P + P + A (25) 

Where in Eqj25| p represents the proton, A represents the interacting nucleus, and p rep- 
resents the antiproton. After the interaction all particles, including the antiproton, continue 
through space until colliding or reacting with other particles. 

3.2 The Tevatron 

Although CERN built the first machine with proton-antiproton collisions [2^ , the highest 
energy model came in the flagship design of the Tevatron at Fermilab. Originally designed 
for fixed target collisions, the main collider ring at Fermilab up into the 70 's could only 
accomplish orbiting energies of about 400 GeV. In 1977 however. Bob Wilson foresaw a 
way of constructing a new ring attached to the existing main ring that would double the 
then current maximum orbiting energy of an injected proton. Essentially he wanted to 
construct a "Super Ring" that would implement the relatively new concept of magnets that 
accomplished superconductivity by supercooling. In creating this Tevatron Super Ring, a 
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proton could be ramped up to an orbiting energy of 1 TeV, then passed through the main 
colhder ring to interact with another proton at 250 GeV [28] . 

After the completion of the Tevatron, initial applications came in ramping a proton to 800 
GeV and smashing it into fixed targets. Soon after CERN showed the possibility of sustaining 
and colliding antiprotons, however, the Tevatron had to switch gears. In 1985, Fermilab 
introduced the Antiproton Source into the collider design that could create antiprotons 
and send them into the Tevatron for collision with protons. All collisions therefore would 
occur and be monitored in the Tevatron itself. Vast improvements have been made in 
antiproton production and peak luminosities since that time. Recently, peak luminosities 
were on the order of 3 x 10'^^ cm~^sec~^. Just in the 5 year span after slip stacking was 
introduced in 2005 [2S], weekly antiproton production went up roughly 8 times (~ 500 x 10^° 
antiprotons/ Week in 2005—)- ~ 4000 x 10^° antiprotons/Week in 2009) [30j. The most current 
data from Fermilab (Figj5]) shows peak luminosity to be on the order of 4 x 10^^ cm~^sec~^ 
[3T| and this is the luminosity we will use to discuss design implementation into our proposed 
system in the next section. 



Tevatron Collider Run I & II Lumincsitv 



If 32 



1131 



s 




1{3|] 



If 29 



^Jt* jtjS\ ,J8l « 



Figure 5: Current Fermilab Luminosity Data |3T 
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3.3 Future Projections 

Although there are many implementations that can be made to improve upon the current 
Tevatron design, we will focus on the a few that could be applied to the implementation of a 
new collider, similar to the Tevatron, to be installed in the VLHC [321 133], tunnel that was 
the focus of Chapter 1. One immediate concern that arises out of the VLHC design we've 
proposed, is the sheer circumference difference between the current Tevatron (6.28 km) and 
the VLHC we've proposed (233 km). If the collider ring circumference was in fact 233 km, 
then collision frequency would go down roughly by a factor of 37 (|^ ~ ^). However, 
because we're talking about a final collision energy of 20 TeV versus the 1 TeV currently 
implemented in the Tevatron, you gain a factor of about a 20 times more energetic collision, 
thus your luminosity would scale roughly as 

Luminosity proposed = {Energy increased x Frequency decreased) x Luminosity current 

X 

Luminosityproposed = (20 x — ) x (4 x 10^^ cm"^sec~^) = 2.16 x 10^^ cm~^sec"-^ (26) 

So it's easy to see that we've basically gone down by a factor of two when simply consid- 
ering putting current designs, with a higher collision energy, into the 233 km collider ring. 
This is a fairly high peak luminosity, but we can still do a bit better by increasing another 
factor, a. Although proton-proton (pp) and proton-antiproton (pp) collision are virtually 
indistinguishable for low mass objects as threshold mass objects are approached, as mass 
and energy increase, there is a verifiable difference in the behavior of the two interactions. 
Simulating three scenarios including: pp at LHC 7x7, upgraded pp at LHC 20x20, and a 
proposed pp at VLHC 20x20, you can see a vast difference in the behavior at higher mass 
collisions in Figj6] |3l]. See also Fig. FigJT] |38j- 

Finally we see that when considering the current LHC parameters applied to pp collisions 
compared to the proposed collisions of pp, at higher mass collisions, cr goes up by a factor 
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Figure 6: Cross Section vs Collision Mass [34] 



of 10. Most importantly, when considering the Event Rate, it must go up by a factor of ten 
as well (a — )■ 5 X 10~^ — )■ 5 x lO^^). This works out great, as the Event Rate goes as. 



Event Rate = {Luminosityproposed x a) 
Event Rate = {Luminosityproposed x lOo") 
10{Event Rate) = ((2.16 x 10^^ cm-^sec"^) x lOa) (27) 

So then by EqJJTJ we clearly have a 10 times higher event rate when considering pp 
interactions at higher mass collisions compared to pp for the same luminosities. So we can 
see that interesting event rates are about the same, but the background collisions are lower. 
Thus a lower cost detector with fewer channels is possible while keeping channel occupancy 
constant. 
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So for high mass objects near theshhold, this colhder, with the current Fermilab p source, 
has double the event rate of the Superconducting Super Colhder (SSC) design. The current 
limitation on the Fermilab p source is the accumulator ring which has a p stacking rate of 
26 X 10^° p/hour. The debuncher ring can supply 40 x 10^° p/hour, so a second p accumulator 
ring used in parallel might help. 
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3.4 Collider Design 

To date there have been many particle colhder built and functioning around the world. 
Many more have been considered, including the Very Large Hadron Collider, or VLHC [32j. 
Previous work at VLHC used roughly 100,000 amps. This current came in the form of a 
magnet construction that entailed a half turn around two bores with a Niobium Titanium 
superconductor at T = 4K with liquid helium in an 8 cm diameter pipe. 

The design could be altered to bring the current consumption down by a factor of 4. 
In the previous PIPETRON design (FigJH]), a centered circular tube of current directs two 
paths for beam collision (FigJH]). By switching to an H-frame dipole construction (Fig JTU]) . 
you would only have one beam path, thus cutting current by a factor of two, to 50,000 amps. 
By implementing a full turn and one bore around the H-frame dipole magnet, only one offset 
current loop would be necessary, thus returning the return current back through the magnet. 
This would buy you another factor of two, bringing current consumption to 25,000 amps. 
With only one offset current loop in the H-frame dipole there could exist the possibility of a 
skew quadrupole field, but nevertheless, a lattice may still be possible [39] . 






Figure 8: Pipetron Schematic. |10] 
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Figure 9: Finished Pipetron Prototype. [1T| H2] 
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Figure 10: H-frame Dipole Schematic |13] (The mitered joints are needed for grain oriented 
sihcon steel, but not with low carbon steel.) 



Another possible design change could come in the form of substituting SuperPower YBCO 
superconductors with T = 25K triple point liquid neon instead of liquid helium [HI HSl |16] . 
By implementing liquid neon and SuperPower YBCO, cooling costs could be minimized. The 
availability of liquid neon is adequate. When considering the heat of vaporization (HOV) 
of each element, neon has an HOV of approximately 1.71 kilo- Joule/mole, whereas helium 
has an HOV of approximately 0.083 kilo- Joule/mole. This means that neon could absorb 
roughly 21 times as much heat per mole as helium in a quench, if the superconductor goes 
normal. Also, considering that both gaseous helium and neon are lighter than air, neon 
might also be usable in a tunnel where oxygen deficiency is a concern. So we should estimate 
the amount of neon needed to accomplish such a collider design. 



V olume collider = 2 X 233 km x 7r(.5 cm)^ = 36.6 million cc (28) 
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Now if we assume that the hquid density of neon is 1.2 g/cc and the gaseous density 
of neon is 0.0009 g/cc, then the gaseous form is roughly 1,300 times less dense than the 
liquid form. Then it is easy to see that collider ring would require 47.6 million liters of neon 
gas. Linde currently sells neon gas for approximately 8 cents/liter, so the approximate cost 
of neon gas would come in at $3.8 million. This cost is pretty good, but the neon would 
necessitate recovery after quenches. Linde currently produces about 300 million liters of 
gaseous neon per year. Their delivery method requires full size tube trailers that hold 3.2 
million liters each, so a project of this size would require a fleet of 15 Kenworth big trucks, 
each hauling a tube trailer to complete the delivery. 

Implementing the SuperPower YBCO Superconductors would be more challenging. Cur- 
rently 0.01 mm X 12 mm 2G HTS ribbon costs $70/meter and has the physical parameters 
to handle 1,500 amps at T = 25K. To get the 25,000 amps needed for our proposed design, 
a minimum of 17 ribbons would be required. Then the total length of YBCO Superconduc- 
tor ribbon would amount to 17 x 2 x 233, 000 m = 8, 000, 000 m. At the proposed cost of 
$70/meter, the total cost of the ribbon would be about $555 million. This cost seems a bit 
pricey and it may turn out that Niobium Titanium and liquid helium at T = 4K are cheaper, 
but the cryogenic load will be far greater at 4K than at 25K. 



26 



4 Muon - Muon 



4.1 Theory 

Muons that are accelerated or stored in collider rings |Tnii8lli9l[50l[5Tll52l[53ll5ll[55ll56l 
1571 [58] will decay. The simplest form of muon decay comes in the form of Eqj29l where a muon 
decays into an electron, an electron-antineutrino, and a muon-neutrino. The alternative 
decay is represented in the second half of Eqj29l in which the antimuon decays into a positron, 
an electron- neutrino, and a muon-antineutrino. These decays happen continuously and will 
eventually contribute an intense disk of neutrino radiation emanating tangentially out of the 
plane of a collider ring [59]. 



/i 1^1, + iye + e 

IJ,^ ^ 17^ + Ue + e+ (29) 

Muons typically have a lifetime of 2.2 /is, however the muons in the scope of this project 
are boosted to energy levels on the order of 20 TeV, so their lifetimes are greatly increased. 
The boost factor 7 is defined in Eq|30]and a sample calculation of the boosted mean lifetime 
is included as well. 



7 



20 TeV 
7 = = 188,680 

' 106 Mev 

Uoosted = X 7 = 2.2 /is X 188, 680 = .415 s (30) 
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A very direct concern comes in the size of the beam that ejects at the earth's surface. 
Essentially this beam width can be approximated by examining the amount of energy each 
neutrino receives from the decay and what amount of that energy is contributed in the 
perpendicular "kick" to the beam path direction. If we approximated each particle in the 
resultant muon decay to receive a third of the 20 TeV muon directed energy, each component 
should travel with about 6.6 TeV. In the case of the neutrino, the perpendicular kick should 
be about 20 Mev, so the appropriate opening angle is illustrated in EqJSTl 



= = 1^ - 3,radian. (31) 

Adjacent 6.6 lev 



4.2 Past Research 

Many reports have explored and approximated the harmful effects of a high energy muon 
collider and several more factors have to be considered and addressed. One such concern 
is the amount of neutrino radiation dosage that exits at the earth's surface. One such 
approximation comes in the form of Eql32l [59] 



DTJASv] = 2.9 X 10-24 X X (^^[TeV])^ 

D[m\ 

If we consider the collider ring to be 300 meters below ground (D = 300 m), the energy 
of the muon to be approximately 20 TeV ( E = 20 TeV), and the constant to be relatively 
accurate, we only have one variable left to define. The number of muons we need per year 
goes down greatly because their lifetimes have gone from roughly 2.2 fis to 0.415 seconds. 
This means that the number of bunches and cycles goes down greatly because fewer bunches 
of muons are needed as the muons created simply live longer. If we consider 2 bunches of 
2 X 10^2 muons per cycle, with each muon having a mean lifetime of 0.4 sees, then we arrive 
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at about 3.1 x 10^° muons per year. Plugging all these parameters into Eql32] will yield 
Eqj3l 



DZtASv] = 2.9 X 10-24 X ^^-^ ^°3oo''J^° ^""^^^ = ^400 [mRem/year] (33) 



This yearly dosage is roughly 24 times more than the federal limit of approximately 100 
mRem/year. 

4.3 Future Projections 

Although there are several approaches to minimizing the yearly dosage of neutrino ra- 
diation, perhaps the most straightforward solution would come in simply spreading out the 
size of the radiated beam. We must first consider a hypothetical circumstance in the present 
approximated beam with an opening angle described in the previous section of this paper. If 
we approximate the beam surfacing at a distance of 250,000 feet away, the beam width that 
surfaced would be approximately 0.75 feet (9 inches). This tightly concentrated beam would 
exceed allowable radiation limits, if considered to interact with a dense surface (concrete 
wall etc.), before showering a bystander, over a long period of time (all day, everyday, for 
years). To spread out this concentrated beam may be a fairly simple fix. 

A roller coaster lattice is implemented in which the FODO (Focused-Off-Defocused-Off) 
quadropoles [60] make the muon bunches oscillate up and down to offset the continuos 
neutrino decay only in a plane tangent to the track. A variation of this design can be 
seen in the helical orbit focusing of protons and antiprotons implemented at the Fermilab 
Tevatron [61]. Essentially, in the helical model you are continuously altering the path of 
the protons and antiprotons in their vertical and horizontal planes to have them spiral past 
one another without ever touching until they reach their specific collision site. This helical 
model is a pretty straightforward addition to an existing collier ring and only requires adding 
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electrostatic separators into the existing collider ring design. 

For muons, path oscillation would only occur in one plane, much like a roller coaster. 
Simply introducing a roller coaster shaped muon stream instead of a planar circular track 
would not have any immediate advantage as the tangental neutrino decay would remain 
relatively the same. The spreading of the otherwise tightly focused neutrino beam would 
come from introducing periodic bumps via short vertical dipoles into the track of the collider 
ring. By adjusting these bumps at random times, you could create a virtually random 
succession of phase shifted beams. This phase shift would ultimately lead to a great spreading 
of neutrino beam interactions as it would spread out, at random, the vertical angle of the 
resultant neutrino decay at all points around the collider ring. If the muons could rise 1 
cm over 20 m, then an angle of 500 /iradians could be achieved. This would be roughly 160 
times higher than the previous 3 /iradians opening angle, more than curing the 24 times 
neutrino radiation dosage problem. To put this in prespective, a dose of 15 [mRem/year] is 
equivalent to eating a banana at every meal. 

4.4 Collider Design 

When considering the construction of a 2 Tesla magnet for a 20 x 20 TeV muon collider 
in a 233 km long collider ring, previous work at VLHC used roughly 100,000 amps. If we 
consider the initial muon injection energy to be approximately 2 TeV and the energy gained 
per revolution of the entire collider ring to be approximately 250 GeV, then the approximate 
number of orbits to reach a overall energy of 20 TeV will be about 72. If we take the muons 
to move near the speed of light, then we see that they make the 72 orbits in about .06 
seconds. This amounts to a 180 degree half cycle if an offset White circuit [62] runs each 
magnet string, then the resultant frequency of the magnets becomes 8.3 Hz. From previous 
research, many superconductors have been used on a small scale at up to 60 Hz [631 E3i ES] . 
Knowing that the eddy current losses go as the square of the frequency, the eddy current 
losses would be about 50 times less at 8.3 Hz as compared to 60 Hz. 
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If we wanted to calculate the eddy current losses within the dipole magnets, such calcula- 
tions are pretty straight forward [661 [57]. Taking ultra low carbon steel 0.5 mm laminations 
to reach 2 Tesla, the dipoles ramp from 0.2 Tesla to 2 Tesla, for a swing of ± 0.9 Tesla. 
Taking the volume of the steel to be 0.25 m x 0.25 m x 233, 000 m = 15, 000 m^ (roughly 
2.5 times as much steel used in the LHC), we can calculate the power consumption of eddy 
current losses. 



Volume x(2x-7TxfxBx 
Power = 



Power 



24p 

15, 000 X (2 X TT X 8.3 x .9 x .0005)^ 



24 X (9.6e-9) 

Power = 36 Megawatts (34) 

But the duty cycle is roughly = 0.3 and not continuous, therefore the 36 Megawatts 
becomes about 11 Megawatts. This consumption is well below the 200 Megawatts being 
discussed for other projects. 

Another power loss must be considered in the form of hysteresis losses. If we consider 
the hysteresis losses in the dipoles which cycle around the BH loop every .4 sees, the energy 
goes as the area of the BH loop divided by An in CGS units. If we approximate the losses 
by a very straight forward method, the Steinmetz Law, we can estimate such losses [HTJ [68] . 

Energy/cycle = (.001 Oersted) x (9000^'^ Gauss) = 2100 ergs/cc (35) 
then plugging Eqj35]in EqJMjwe get. 



p^^^^ _ (15,000 m^) X (10^ cc/m^) x (2100 ergs/cc) x (10"^ ergs/joule) 

.4 sees 

Power = 8 Megawatts (36) 
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Again we're still well below a 200 Megawatt site power limit. It also should be noted 
that substituting grain oriented silicon steel for very low carbon steel would lead to a large 
decrease in both eddy current and hysteresis losses. But, in exchange, the maximum field 
would decrease from 2 Tesla to 1.8 Tesla, therefore leading to a 10% lower beam energy. But 
losses in annealed 0.005% ultra low carbon steel [691 EO] appear to be low enough. The muon 
collider magnets ramp quickly, so eddy current and hysteresis losses are a concern. The pp 
magnets ramp slowly. One may want to ramp the e~^e~ magnets quickly to go from 20 to 
44GeV in two orbits. This minimizes the time for the transverse mode coupling instability 
to cause trouble. The RF to do this is already in place and will be used at higher energies 
to make up for synchrotron radiation losses. 

Now we can also see how the luminosity, L, scales with energy and collision frequency. 
L = fnNiN2/A where / is the revolution frequency, n is the number of bunches in one beam, 
Ni is the number of particles in each bunch, and A is the cross sectional area of the beam 
at the interaction point. Also it is important to note that A = AnaxCy, where cx^ and Cy 
are transverse sizes. The transverse angular beam spread decreases linearly with energy as 
seen in o" = a/cat/?*/ /S^y where e is the normalized emittance in mm-mrad, /3* is the collision 
region focal length, /3 = v/c, and 7 = 1/^/1-/3'^. So for constant focal length, /?*, the 
decrease in collision frequency in a larger ring is exactly canceled by the higher muon energy, 
which goes as 7, and the luminosity stays the same. So at least then the luminosity doesn't 
go down. 

We must also consider the overall energy of the beam once it reaches the maximum muon 
energy of 20 TeV. This calculation is also indicative of what the overall beam power would for 
the specific muon lifetime of ~0.4 sees. This specific lifetime essentially means that the muon 
injector must replenish the muon supply every 0.4 sees, giving it an operating frequency of 
approximately 2.5 Hz. Then, 
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Joules 

Total Energy max = 2 x [NMuons) x {Emax) x — — 

eV 

Total Energymax = 2 x (2e^^) x (20e^^) x 1.6e"^^ = 1.28 x 10^ Joules. (37) 



Then Power goes as, 



Power = Emax x / 

Power = 1.28 X 10^ Joules x 2.5 Hz 

Power = 32 Megawatts. (38) 

So it's easy to see that maximum power is still below the 200 Megawatt site limit, even 
when efficiency is added in. The last calculation to make for the proposed collider design 
will be a model that replicates the muon survival rate when ramping up to the maximum 
energy of 20 TeV. A few parameters must be defined, but essentially it acts as a product 
series of an exponential function. The muon survival rate can be defined as 

Orbits 

SURVIVAL = Y[ exp 

Af=l 

So for the collider design proposed, we want to ramp up to 20 TeV, in 72 orbits. We 
set the Einjection to 2 TeV. We also set the energy gained per orbit around the collider ring 
of the muon, Eorbit, to be .25 TeV. Now considering we are talking about the proposed 233 
km collider ring construction, the respective radius, R, would be 37,100 m. Now we have 
all the parameters we need to evaluate the product summation of the series, but evaluating 
the series itself is not such a trivial calculation. For this reason, I created a Mathematica 
program to evaluate this product. The program is shown in FiglTTl 

Therefore evaluating aval = 37100 for the ring radius in meters, bval = 2000 for the 
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injection 



+ {EorMt N)] CT 



Survival%. 



(39) 
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Einjection in GeV, and the cval — 250 for the Egrbit in GeV, we calculated the survival rate 
to be approximately 



72 

SURVIVAL = Yl 

N=l 
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Figure 11: Mathematica NoteBook Example. 
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71.58%. 



(40) 
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5 Impacts on Society 



5.1 Future Projections 

Robotic tunneling may be advanced, particularly in the areas of rock removal |71]. and 
shotcrete tunnel lining. Advancement of the use of high temperature superconductors at 25 
Kelvin may lead to less expensive magnetic resonance imaging machines in medicine. It is 
much easier to reach a temperature of 25 Kelvin than 4 Kelvin. The multi-megawatt proton 
accelerators such as fixed field, alternating gradient FFAG rings needed to produce muons 
may also be able to fission thorium [72j to produce copious amounts of electricity. Thorium 
reactor fission stops automatically in the absence of the proton beam. 
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Conclusion 



Ten years ago in 2001, the cost of a 233 km circumference tunnel in Northern lUinois was 
estimated to cost $2.8 biUion including 30% contingency. Since then inflation has increased 
prices but more automation has been added to tunneling in pulling tunnel boring machines 
forward and in placing rock stabilization bolts. This leads me to beheve that construction 
costs would be directly on point with figures proposed over a decade ago, and perhaps even 
a httle lower. 

A 500 GeV center of mass e^e~ collider with a luminosity of 10^^ cm~^ s~^ appears to be 
possible. A crab waist crossing as developed for the next generation B factories is employed 
and shows great promise for high luminosity by exploiting low emittance. Low emittance 
bunches from precision damping rings for the proposed International Linear Collider (ILC) 
are used as well the short vertical focal length ILC collision region optics. Forty times 
less installed RF acceleration is required in the ring as compared to a linear collider. Low 
coercivity grain oriented silicon steel is used for dipole magnets. 

A 40 TeV pp fits in the tunnel with 2 Tesla H-frame iron dipoles. For many high-mass 
objects, the pp cross section is an order of magnitude larger than the pp cross section. Single 
bore iron dipoles permit high temperature superconductor ribbons running at liquid neon 
temperatures. Advances have been made in p production at Fermilab. A detector with more 
physics and fewer background interactions per crossing requires fewer channels and therefore 
as a consequence costs less. 

A n'^ ii~ collider approaching a center of mass energy of roughly 40 TeV fits in a 233 
km circumference tunnel. The same dipoles as used for pp can be employed. Eddy current 
and hysteresis losses are shown to be reasonable. The ramping rate is low enough so that 
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superconductor can be employed. With 250 GeV of RF the muons are accelerated in 72 
orbits. Muon survival ~ 72%. A phase shifting, FODO, roller coaster lattice similar to the 
Tevatron helical lattice is used to mitigate the neutrino radiation problem down to a level 
equivalent to eating one banana per meal. The vertical spread of the neutrinos is simply 
increased. Half the RF of the ILC is used to produce a lepton collider with 80 times the center 
of mass energy of the ILC, if muons can be cooled [SH] at a reasonable cost [75| Beam 
power for the muon and electron machines [16] is similar, 32 Megawatts vs. 21.6 Megawatts, 
respectively. 
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